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Surface bound states in a discrete-lattice model of a d x i_ y i cuprate superconductor are shown 
to be, in general, coherent superpositions of an incoming excitation and more than one outgoing 
excitation, and a simple graphical construction based on a surface Brillouin zone is developed to 
describe their nature. In addition, a momentum-dependent lifetime contribution to the width of 
these bound states as observed in tunneling experiments is derived and elucidated in physical terms. 

PACS numbers: 74.20.-z, 74.25.Jb, 74.50. +r, 74.80.fp 
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One of the consequences of the change of sign of the 
gap function in the d x 2_ y 2 model of the high-temperature 
cuprate superconductors is the existence of zero-energy 
surface bound states which are believed to be responsible 
for the zero-bias conductance peaks (ZBCP's) ob- 
served fj|-|l4| in the in-plane conductance of tunnel junc- 
tions. Further support for the idea that surface bound 
states are responsible for the ZBCP's comes from the 
study of the splitting of the ZBCP's in a magnetic field 
[ j||n] . Furthermore, the splitting of the ZBCP's observed 
[ ll| at low temperatures in zero magnetic field provides 
evidence Hwj] for the existence of the so-called d + is sur- 
face state [ptpjjlq] (where the order parameter acquires an 
s-wave component at certain surfaces which is 7r/2 out of 
phase with a suppressed d-wave component). 

In current models [jl]-[7|, the surface bound state is a 
coherent combination of an incoming excitation and an 
outgoing excitation whose existence is conditional on the 
gaps of the two excitations having opposite signs. Since, 
in a d x 2_ y 2 superconductor, the energy gap has changes 
of sign imposed by symmetry, surface bound states are a 
natural consequence of d x 2_ y 2 symmetry, which is one of 
the reasons for their current interest. 

Here we study a discrete-lattice model which gives sur- 
face bound states with new features. In particular there 
are surface bound states which are coherent combinations 
of an incoming excitation with more than one outgoing 
excitation provided all outgoing excitations have gaps of 
the same sign, which must be different from the sign of 
the gap of the incoming excitation. The existence of sev- 
eral conditions (one for each outgoing excitation) reduces 
the number of surface bound states, as can be seen below 
in detailed results for (210) and (110) surfaces. An im- 
portant aspect of the work below is the introduction of a 
simple graphical construction, based on the application 
of Bloch's theorem to translations parallel to the surface 
and the idea of a surface Brillouin zone, which elucidates 
the more complex structure of the surface bound states 
in a discrete lattice model. 

We also develop a theoretical expression for the ZBCP 



contribution to the tunneling conductance, and in partic- 
ular elucidate a mechanism of a momentum-dependent 
broadening of the surface bound states which must in- 
evitably accompany attempts to observe these states in 
tunneling experiments. This lifetime broadening comes 
about because a quasiparticle in a surface bound state on 
the superconducting side of a normal to insulating to su- 
perconductor (NIS) tunnel junction will ultimately tun- 
nel into the normal metal. If the probability of tunneling 
through the insulating barrier is P (ky) (which depends 
on the momentum parallel to the surface, ky) then the 
lifetime r(k||) will be given by t" 1 (kjj ) = P(ky)/(ky) 
where /(ky) is the so-called attempt frequency, i.e. the 
frequency with which the quasiparticle returns to the sur- 
face. (This is a familiar result from textbook studies 
of the alpha decay of nuclei.) To calculate the attempt 
frequency, recall from previous studies of surface 

bound states that the distance the (exponentially decay- 
ing) bound state wave function extends into the super- 
conductor is £(ky) = ftt>j_(k||)/A(kn) where A is the gap 
and w_l is the quasiparticle velocity normal to the surface; 
the attempt frequency is thus / ~ v_l/(2£) = A/(2h). 
Thus, the energy width of the surface bound state is 
T - h/r, i.e. 



r(ky) = iP(ky)A(ky). 



(1) 



This result has a significant dependence on ky and is the 
width which will be relevant for sufficiently pure super- 
conductors having sufficiently flat surfaces. 

Other articles have taken into account the broadening 
of the surface bound states through impurity scattering 
by adding a momentum-independent imaginary part to 
the energy |]Jl^], and the broadening effect of surface 
roughness has also been considered ||. Also, the lifetime 
effect elucidated here is presumably implicitly present in 
the numerical plots of tunneling conductance presented 
by some other authors but its physical origin as 

described above does not appear to have been commented 
upon, nor have explicit expressions for this lifetime and 
its momentum dependence been given. 
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For an NIN junction, we find for the conductance per 
surface unit cell, Gjv, the result 



G N = 2(e 2 /h)(P(k«)} 



(2) 



where the angular brackets indicate an average over the 
wave vectors k|| . This has the familiar and intuitively 
appealing form of a Landauer formula flq| . Clearly a 
measurement of Gn provides information relevant to an 
estimate of the width r(ky) through Eq. (0). 

In the case of an NIS junction, we show below that 
the conductance per surface unit cell associated with the 
surface bound states, and which gives rise to a zero bias 
conductance peak (ZBCP), is 



Gzbcp = 4— 
h 



r 2 (k,i) 



(e^-£: B (k||)) 2 + r 2 (k|| 



(3) 



where r(ky) is given by Eq. (0) , and E B (ky ) is the energy 
of the bound state. The term in the angular brackets is 
the probability of Andreev reflection (clearly a resonant 
process) of an electron incident on the junction from the 
normal side [and having energy less than A(k||)]. Inter- 
estingly, when the voltage satisfies eV — Eb(\s-\\), the 
probability of Andreev reflection for an electron at that 
value of k|| is unity. The factor of 4 in this equation 
comes from one factor of 2 accounting for both electron 
and hole processes, and a second factor of 2 coming from 
the fact that the Andreev reflection of a single electron 
results in two electrons being transferred across the bar- 
rier into the superconductor. Also, since £ , s(k||) is an 
odd function of k||, Gzbcp is an even function of the 
voltage. While Eq. (||) is limited in that it ignores the 
conductance associated with the usual propagating exci- 
tations in the superconductor, it is useful in that it gives 
a simple, intuitive and detailed expression for the ZBCP, 
which is a separate and clearly identifiable part of exper- 
imental conductance versus voltage data. Also, although 
Eq. (||) is derived below in detail only for the case of a 
(110) surface, results of the same form are obtained for 
different models of the surface with different results for 
Eb, and it is therefore tempting to speculate that the 
result is of general validity. 

We now outline some of the technical details, beginning 
with the description of the possible surface bound states 
for the (210) surface shown in Fig. [l](a). We describe the 
superconducting state by the discrete-lattice Bogoliubov- 
de Gennes equations 



E 



-t(x,x') A(x,x') 
A*(x,x') £(x,x') 



U{x') = EU(x) 



(4) 



where U(x) is the two-component wave function and x la- 
bels the ion positions. For simplicity we consider nearest- 
neighbor hopping only [i(x, x + a) = i(x, x + b) = 
t] and a e?-wave order parameter [A(x, x + a) = 
-A(x, x + b) = A]. 



Note that the crystal, including its surface, is invari- 
ant with respect to a translation parallel to the surface 
through a distance \/5a. Thus, from Bloch's theorem, 
the eigenstates of Eq. (0) can be written in the form 



(5) 



where \k y \ < n / \-\fba), and m and n are integer ionic row 
and column indices as indicated in Fig. [l]. Substituting 
Eq. (§) into Eq. (@) yields a one-dimensional equation 



T k y (n, n')U k y {n) = EU k y (r 



(6) 



for each k y . Note that the nearest-neighbor interaction 
between the ions in the two dimensional lattice couples 
the ionic column n to columns n! = n — 2, ?i — 1, rt, n + 1 
and n + 2. Thus Eq. (|) is a fourth order difference equa- 
tion which will have four linearly independent solutions. 

Before describing the solutions of Eq. (0) it is neces- 
sary to introduce the idea of a surface Brillouin zone. 
The normal Brillouin zone for the square lattice of Fig. 
[j] is the square in Fig. || containing the set of wavevec- 
tors |fc |, \kb\ < 7r I a . The surface Brillouin zone is the 
rectangle in Fig. ||(a) containing the set of wave vectors 
\k x \ < Vbir/a, \k y \ < n/(\^5a). These two Brillouin zones 
are equivalent in the sense that a complete set of bulk 
states can be equivalently described as having wave vec- 
tors in either Brillouin zone. However, from Eq. (0), an 
arbitrary surface eigenstate is most simply described as 
a linear combination of bulk states whose wave vectors 
lie on a line of fixed k y in the surface Brillouin zone. 
The Fermi surfaces in Fig. are drawn not as calculated 
from our model, but as determined from angular resolved 
photoemission studies |l7[]l8f| from YBa2Cu3 06+a; and 
Bi 2 Sr 2 CaCu208+a:- 

The bulk solutions that can participate in the forma- 
tion of a surface bound state with a given k y have energies 
less than the minimum energy gap A(k) where k is one 
of the wave vectors where a line of fixed k y cuts the Fermi 
surface. As in previous work Jl]-0] such solutions can be 
found as an expansion in powers of A/e^, tp being the 
Fermi energy, and, to lowest order in A, are 



U- y (n) 



E ± ifl,- 



e mfc ix a/v'5 e - 



(7) 



Here ki = (k y ,ki X ) must be on the Fermi surface; also 
Ai = A(ki), Qi = y/A'j-E 2 , and = / {Vhh\v ix \) 
with Vi X being the x component of the normal-state elec- 
tron velocity. The upper and lower signs in Eq. (Q) cor- 
respond to Vi x < and Vi X > 0, respectively. From 
Fig. ||(a) it is seen that for many values of k y , a line of 
constant k y intersects the Fermi surface at four distinct 
points, giving four values ki X , i= 1 ... 4, and hence the re- 
quired four linearly independent solutions. In these cases 
the surface bound state solutions will have the form 
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U k v(n) = Z i C i U- v (n), (8) 

where only one incoming wave (i.e. V\ x < 0) is included, 
for which C\ is taken to be unity, as are two outgoing 
waves for which Ci and C3 (in general of comparable 
magnitude) are to be determined. The boundary con- 
ditions determining C'2, C3 and the surface bound state 
energy E are Eqs. (^|) for n = and n — 1, which yield 

U k *(-1) = U k y(-2) = 0. (9) 

A surface bound state with E = is found if the outgo- 
ing wave gaps A2 and A3 have the same sign, and are 
opposite in sign to the incoming wave gap Ai. 

If a line of constant k y intersects the Fermi surface at 
only two points in the surface Brillouin zone, only two 
solutions (say with gaps Ai and A2) are of the form of 
Eq. (Q), and a further two solutions must be sought nu- 
merically. As in Eq. (]?]), these solutions are proportional 
to exp(iqn) where q is complex. In Eq. ((7J), however, 
Im(q) <C 1 by virtue of the approximation (Ai/ep) <C 1, 
whereas in the present case the two additional solutions 
will have Im(q) ~ 1 and are thus localized close to 
the surface. All four solutions are included in Eq. (||), 
which gives a surface bound state with E = provided 
AiA 2 < 0. 

Applying the conditions just established shows that 
the only surface bound states which can occur on a (210) 
surface are those made up from incoming waves whose 
wave vectors fall on the Fermi between the two dashed 
lines in Fig. ^. There are relatively few of these bound 
states because having to satisfy a condition on each of two 
outgoing waves is more restrictive than having to satisfy 
a condition on only one outgoing wave, as in other studies 




We have also determined the surface bound states 
for a (110) surface (as in Fig. |l|(b) but including only 
ion columns n > 1), and with the model extended 
to include a second-nearest-neighbor hopping interac- 
tion i(x, x + a + b) plus related terms derived from the 
tetragonal symmetry (which gives a more realistic Fermi 
surface, comparable to that shown in Fig. ||). The formal 
aspects of the problem, including the boundary condition 
and the conditions for the existence of surface bound 
states, are analogous to the problem studied above of 
a (210) surface for the case of nearest-neighbor interac- 
tions. Applying these conditions leads to the conclusion 
that surface bound states with zero energy exist only for 
those values of k y which lie between the two dashed lines 
in the surface Brillouin zone in Fig. ||(b). 

We now describe the calculation of the ZBCP tunnel- 
ing conductance for the NIS junction at the (110) surface 
shown in Fig. 0(b), restricting our discussion to the fea- 
tures particular to our model, since the general features 
of such calculations are well established |,|,|. The 10ns 
of the column n = in Fig. pVb) have an additional 



repulsive potential v > and thus represent an insu- 
lating barrier. We return, for simplicity, to the case of 
a nearest-neighbor only hopping interaction, t (assumed 
to be the same for all pairs, metal, insulator, or su- 
perconductor) giving a normal-state quasiparticle energy 
e(k) = -Atcos(q x )cos(q y ) where q x , v = k x , y a/V2. 

Assuming for the moment that the metals on both sides 
of the insulating barrier are normal, the probability of 
transmission of an electron on the Fermi surface across 
the barrier is easily found to be 

P(k y ) = l6(t/v) 2 (cos 2 q y - v 2 ). (10) 

The chemical potential has been written in terms of v 
by fi = Atv (so that v = and v — 1 correspond to 
half-filled and full bands, respectively) and q x is deter- 
mined in terms of q y for electrons on the Fermi surface by 
cosq x = —v I cosq y . As expected, it is easiest for electrons 
incident normally on the barrier to be transmitted, and 
the probability of transmission goes to zero for electrons 
travelling parallel to the barrier. 

As described above, a constant nearest-neighbor gap is 
assumed for ions in the superconductor except at the sur- 
face where we add an additional gap (5A(x, x + a) = A a 
and 5A(x, x + b) = Af,; this additional gap is nonzero 
only for x lying on the superconducting surface, i.e on 
the column n = 1. The definitions of A s andA^ (both 
real) by A a + Ab — 2iA s and A a — A& = 2Ad corre- 
sponds to the existence of a d + is state at the surface, 
as suggested in ||p|,P,|7|,|l5|] . Also, although our calcu- 
lation is not self-consistent, the expected suppression of 
the d-wave component of the gap at the surface can be 
mimicked by an appropriate choice of A^. 

Making use of the experimental fact that the width of 
the ZBCP is much less than A, and also taking £« A, 
we find that the tunneling conductance is given by Eq. 
(||), where the probability of Andreev reflection of an 
electron incident on the junction from the normal side is 
given by the expression in angular brackets in this equa- 
tion. Also, T(k y ) is given by Eq. (|l|) with P(k y ) given 
by Eq. (0) and A(k y ) = AA(cos 2 q y - v 2 ) x l 2 tanq y . The 
energy of the surface bound state, which is independent 
of Ad, i.e. of the d-wave gap suppression at the surface, 
is given by 

E B = -32z/A(A a /t)[l - (v 2 /cos 2 q y )]tanq y . (11) 

Fig. H shows the tunneling conductance as evaluated 
from Eq. ([}]) for several values of the parameters w — 
l6A(t/v) 2 and s — 16A(A s /i) giving the magnitudes of 
the bound state width and shift, respectively. (Note that 
so long as w, s and eV have the same units, these units 
are irrelevant.) In comparing the two curves for which 
there is no surface s + id state, i.e. s = 0, note that the 
momentum dependence of T in the w = 1 case gives a 
sharper cusp-like behavior at the origin relative to the 
case r = 0.5 and independent of momentum. There is 
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also a weak cusp-like behavior at the origin in the broad 
s = w = 1 curve, reflecting the fact that the low energy 
bound states have a low attempt frequency and are thus 
relatively narrow in energy. (Note that zero energy peaks 
in a broad ZBCP have been observed [|llj.) Also, when 
w <C s, the conductance gradually increases when V is 
increased from zero, and then drops rapidly to zero for 
eV greater than the maximum bound state energy. 

In conclusion, we note that the analysis of a discrete- 
lattice model of surface bound states has uncovered new 
behavior, namely bound states made up of coherent su- 
perpositions of multiple excitations, and that a simple 
graphical approach based on the idea of a surface Bril- 
louin zone helps to sort out the detailed structure of the 
bound states. Furthermore a qualitative "alpha decay" 
type of picture (backed up by detailed calculation) is 
given of the momentum-dependent width of the surface 
bound states which would be observed in sufficiently pure 
superconductors with sufficiently flat surfaces. Finally, 
a formula [Eq. (||)] for the ZBCP conductance is given 
which has the possibility to be general in applicability. 

Helpful discussions with M.E. Zhitomirsky are ac- 
knowledged, as is the support of the Natural Sciences 
and Engineering Research Council of Canada. 



FIG. 1. A (210) surface (a) and a (110) NIS tunnel junction 
(b). In (b) the solid (open) circles, and the solid squares rep- 
resent ions in the superconductor (insulator), and the normal 
metal, respectively. The ionic columns and rows are labelled 
by n and by m. 

FIG. 2. The surface Brillouin zone for a (210) surface (a), 
and a (110) surface (b). The thick- and thin-line portions of 
the Fermi surface (sections of rounded squares) correspond to 
positive and negative gap functions, respectively. 

FIG. 3. The ZBCP for a (110) NIS junction as a function 
of the shift s and the width w of the surface bound state. 
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